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We study scalar, electromagnetic and gravitational perturbations of a Reissner-Nordstrom-anti- 
de Sitter (RN-AdS) spacetime, and compute its quasinormal modes (QNM's). We confirm and 
extend results previously found for Schwarzschild-anti-de Sitter (S-AdS) black holes. For "large" 
black holes, whose horizon is much larger than the AdS radius, different classes of perturbations 
are almost exactly isospectral; this isospectrality is broken when the black hole's horizon radius is 
comparable to the AdS radius. We provide very accurate fitting formulas for the QNM's, which are 
valid for large black holes having charge Q < Q ex t/3 (Qext being the extremal value of the charge). 
Electromagnetic and axial perturbations of large black holes are characterized by the existence 
'^f^ ' of pure-imaginary (purely damped) modes. The damping of these modes tends to infinity as the 

| black hole charge approaches the extremal value; if the corresponding mode amplitude does not tend 

, to zero in the same limit, this implies that extremally charged RN-AdS black holes are marginally 

^s^j ■ unstable. This result is relevant in view of the AdS/CFT conjecture, since, according to it, the AdS 

QNM's give the timescales for approach to equilibrium in the corresponding conformal field theory. 
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INTRODUCTION 



(N 

Quasinormal modes (QNM's) are well known to play an important role in black hole physics. They determine the 
late-time evolution of fields in the black hole exterior, and numerical simulations of stellar collapse and black hole 
collisions have shown that in the final stage of such processes ( "ringdown" ) the quasinormal modes eventually dominate 
y—( . the black hole response to any kind of perturbation. For these reasons, QNM's of black holes in asymptotically flat 
spacetime have been extensively studied for more than thirty years (for comprehensive reviews see [HQ))- 

Perturbations of non- asymptotically flat spacetimes have also been studied for years in a cosmological context. For 
example, perturbations of Schwarzschild-de Sitter spacetimes were considered in 0] , and the corresponding QNM's 
were computed in p|. An analogous study for Reissner-Nordstrom-de Sitter perturbations was performed in p|. 

In the last few years, the anti-de Sitter (AdS) conformal field theory (CFT) conjecture has led to an intensive 
5^ ' investigation of black hole QNM's in asymptotically AdS spacetimes. According to the AdS/CFT correspondence 
a large static black hole in asymptotically AdS spacetime corresponds to an (approximately) thermal state in 
the CFT. So the time scale for the decay of the black hole perturbation, which is given by the imaginary part of its 
QNM's, corresponds to the timescale to reach thermal equilibrium in the strongly coupled CFT 0. The computation 
of these time scales in the CFT is a difficult task. Therefore, black hole perturbation theory has turned out to be, 
quite unexpectedly, a useful tool to compute CFT thermalization time scales. 

Quasinormal modes in AdS spacetime were first computed for a conformally invariant scalar field, whose asymptotic 
behaviour is similar to flat spacetime, in Subsequently, motivated by the AdS-CFT correspondence, Horowitz and 
Hubeny made a systematic computation of QNM's for scalar perturbations of Schwarzschild-AdS (S-AdS) spacetimes 
Their work was extended to gravitational and electromagnetic perturbations of S-AdS black holes in ^(j- A 
reconsideration of the correct boundary conditions for the QNM problem was presented by Moss and Norman 
who considered gravitational perturbations of both Schwarzschild-de Sitter and S-AdS spacetimes. Finally, the study 
of scalar perturbations was extended to the case of Reissner- Nordstrom- AdS (RN-AdS) black holes in [l^ . 

Frequency-domain studies of de Sitter and anti-de Sitter black hole spacetimes have lately been complemented by 
time evolutions. In |13| and |14| . the authors used a time evolution approach to study tails in Schwarzschild-de Sitter 
backgrounds from a cosmological perspective. Motivated by the AdS/CFT correspondence, analogous time evolutions 
have recently been performed for RN-AdS backgrounds in [l5| . 

Numerical methods used so far to compute QNM's in AdS backgrounds tend to break down, in the S-AdS case, as 
the black hole horizon becomes much smaller than the AdS radius 0. Therefore, special attention has been devoted 
to the small-black hole limit, using both a frequency domain [T^ | and a time domain approach |17| . 

Our aim in this paper is to present a comprehensive study of scalar, electromagnetic and gravitational perturbations 
of RN-AdS black holes in the frequency domain. We verify (or sometimes disprove) and extend many results obtained 
in the aforementioned papers. We also compute for the first time electromagnetic and gravitational QNM's for 
RN-AdS black holes, building on work done previously in the cosmological context for Reissner-Nordstrom-de Sitter 
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spacetimes 

The plan of the paper is as follows. In section [21 we introduce our notation for the background metric, display 
explicitly the wave equations describing scalar, electromagnetic and gravitational perturbations, and briefly present 
our numerical methods and conventions. Details on the numerical methods are deferred to the Appendix. In section 
13 we first summarize and extend some well-known results for perturbations of S-AdS black holes, considering all three 
cases (scalar, electromagnetic and gravitational perturbations); then we use the uncharged case as a a starting point 
to discuss perturbations of RN-AdS black holes, and present most of our new numerical results. The conclusions and 
a discussion follow. 



2. PERTURBATION EQUATIONS IN A RN-ADS BACKGROUND 

2.1. Background 

We are interested in scalar, electromagnetic and gravitational perturbations of a Reissner-Nordstrom anti-de Sitter 
metric, 

dv 2 

ds 2 = f(r)dt 2 - jj^r - r 2 (d9 2 + sin 2 6d(f> 2 ), (1) 

where f(r) = A/r 2 and 

r 4 

A = r 2 -2Mr + Q 2 + —r. (2) 
R z 

In the previous formulas M is the black hole mass, Q the black hole charge, and R is the AdS radius. The black hole 
mass is related to its charge Q and horizon radius r + by the relation 

while the black hole's Hawking temperature is given by 



T _ l-Q 2 /r 2 + + ir 2 + /R 2 _ 

47T?'+ 



It will be useful in the following to introduce a tortoise coordinate r* , defined in the usual way by the relation 

dr A 



(5) 



The extremal value of the black hole charge, Q ex t, is given by the following function of the hole's horizon radius: 

QL = r\ (1 + ir 2 + /R 2 ) . (6) 



2.2. Scalar perturbations 



The scalar perturbations of the RN-AdS metric have been studied in so we only give the final result. Separating 
the angular dependence of the perturbations and introducing a differential operator 

the radial part of the perturbations obeys a wave equation 

A 2 Z = VZ 

with potential 



V = f(r) 



A 



2 T 



(7) 

(8) 
(9) 



R 2 r 2 



r+r 2 



2Q 2 



2r 



Here / is the usual spherical harmonic index. In the limit Q = this potential reduces to the four-dimensional limit 
of the scalar potential considered by Horowitz and Hubeny 9] . 
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2.3. Electromagnetic and gravitational perturbations 



The electromagnetic and gravitational perturbation equations for Reissner-Norstrom-de Sitter black holes were 
considered in 5] following the track outlined in [2(j ■ Once again, their derivation is carried over with no difficulty to 
the RN-AdS spacetime, so we only give the final expressions for the perturbation equations. After separation of the 
angular dependence, the axial perturbations (which we will denote by a superscript "-" ) are described by a couple of 
wave equations: 



v~z: 



Defining 



and 



(/-!)(? + 2) 



Pi = 3M+ (9M 2 + 8nQ 2 ) 1/2 , 
p 2 = 3M - (9M 2 + 8nQ 2 ) 1/2 , 



the relevant potentials are 



vr 



A 



2(n + l)r - Pi (l + ^)] = 1, 2; j). 



(10) 

(11) 
(12) 

(13) 



Similarly, the radial parts of the polar perturbations (that we will denote by a superscript "+") obey the couple of 
equations: 



K 2 Z+ = V+Z+ 



where 



V 2 + 



A 
A 



U+-(pi-P2)W 

u- l -{ Pl - P2 )w 



(14) 
(15) 



Here U and W are given by 



A It 2r 6 

W = 5 (2nr + 3M) + — [ nr + M - — 7 

rw w \ R- 

( 2r 3 \ 2nA 

U = (2nr + 3M) W + [w - nr - M + — ^ , 

V R I w 



(16) 



and the function w is defined by the relation 



w — nr + 3M 



2Q 2 



(17) 



In the limit Q = 0, the potentials V 1 and reduce to the corresponding potential for pure electromagnetic 
perturbations of Schwarzschild-AdS black holes, equation (5) in fioj ]: 



Vem = A 



1(1 + 1) 



(18) 



Furthermore, in the same limit, the potential V 2 reduces to the potential for pure axial gravitational perturbations 
of Schwarzschild, equation (13) in [Tol ]: 



A 

V dd — o - 



1(1+1) 



6M 



(19) 



and the V 2 + potential reduces to the potential describing pure polar gravitational perturbations of Schwarzschild, 
equation (17) in [Toj : 



2A\ 9M 3 + 3n 2 Afr 2 + n 2 (l + n)r 3 + 3Af 2 (3nr + 3r 3 /R 2 ) 



(3M + nr) 2 



(20) 
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2.4. Numerical method 

In AdS backgrounds, QNM's are denned as solutions to the relevant wave equations characterized by purely ingoing 
waves at the event horizon and vanishing perturbation at radial infinity. We used two different numerical methods to 
solve this eigenvalue problem. The first has been discussed at length in previous papers @,E3>0j an( A we briefly recall 
it using our notation (and covering the slightly more general case considered here) in Appendix^] As a further check, 
we also extended to the axial perturbations of RN-AdS black holes the Frobenius method used in ^lj f° r S-AdS black 
holes, developing on previous work for black holes in asymptotically flat spacetime [Tftl IT^|. This extension (which 
essentially involves the use of a four-term recurrence relation instead of a three-term recurrence relation) is presented 
in Appendix [2| 

The two methods yield numerical results which agree to machine accuracy. In both formulations of the eigenvalue 
problem, we numerically search for roots in the plane (lvr, lot) using Muller's method. For pure-imaginary roots, it is 
sometimes convenient to resort to Ridder's method (for a description of both root-searching techniques, see e.g. |2l|). 
The Frobenius method is computationally much faster, and it is also easier to implement. However, at variance with 
the asympotically flat case, axial and polar perturbations in anti-de Sitter backgrounds are not isospectral [Toj . A 
similar isospectrality breaking was found for perturbations of charged, dilaton black holes in .22]. As a consequence, 
polar perturbations, which are described by the rather messy potentials (|15|) . cannot simply be related to axial 
perturbations through a Chandrasekhar transformation of the radial perturbation variables |2C| . So we only used the 
Frobenius method as a check when computing modes associated to the axial potentials Ijl^jl . 

From now on, in presenting the numerical results, we will set the AdS radius R = 1. It is also useful to normalize 
the charge to its extremal value, defining 

Q - Q/Qext- (21) 

Finally, we set to — lor — iluj, in order to have positive values for the imaginary parts of QNM frequencies. In the 
next sections we will first recall the properties of S-AdS QNM's, verifying and extending results found in [ij. Ho|. 
Then we will present our results for charged black holes, discussing both the scalar modes studied in ^1 an( i the 
electromagnetic and gravitational perturbations, which, to our knowledge, are studied here for the first time. In both 
cases we will first discuss the "ordinary" QNM's, and then the purely damped QNM's (whose real part lor = 0), 
which show rather peculiar features deserving special attention. 

3. NUMERICAL RESULTS 
3.1. Schwarzschild-anti-de Sitter black holes 

3.1.1. "Ordinary" quasi-normal modes 

Consider a sequence of uncharged black holes in anti-de Sitter space, and fix units such that the AdS radius R = 1. 
Then QNM frequencies only depend on the holes' horizon radius r + , and can be plotted as curves in the complex 
(ljr, lji) plane. Such curves are shown for the fundamental scalar mode with I = (continuous line), for the first 
non-purely damped axial mode with I = 2 (dashed line) and for the fundamental polar mode with I = 2 (dotted line) 
in figure n Comparing the dashed and dotted lines in the figure, it is apparent that, as mentioned in the previous 
paragraph, in a S-AdS background the isospectrality between even and odd perturbations (which is a property of 
asymptotically flat black hole spacetime perturbations) is broken. However, it is approximately restored in the large 
lo and in the small black hole limit (for a discussion of this property, see section III C in ^3). Scalar perturbations are 
isospectral with odd and even parity gravitational perturbations in the large black hole (large lo) limit. As yet, there 
is no analytical proof of this rather surprising feature, which we have numerically verified both for the fundamental 
mode and for the overtones. 

In general, our results are in perfect numerical agreement with those presented in jjj ITol ] . In their treatment of 
scalar perturbations, Horowitz and Hubeny found that, for large black holes and I = 0, the real and imaginary part of 
QNM frequencies scale linearly with the temperature. In four dimensions, they predicted that the fundamental mode 
frequencies are well fitted by the relations lor = 7.75T, loj = 11.16T. They also found that, in general, the overtones 
are equally spaced, obeying relations of the form lor{ti) ~ 54+ 131n, u>i(n) ~ 41 + 225n, where n is a positive integer, 
not to be confused with the coefficient n defined in equation 111|) . As / increases, loi decreases and lor increases, but 
the /-dependence of the modes is extremely weak. 

Our numerical results confirm all of these statements, and show that they apply also to axial and polar gravitational 
perturbations. Table[I]shows the fit parameters for the linear temperature-dependence of the lowest QNM frequencies 
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TABLE I: Linear scaling of the real and imaginary part of QNM frequencies with the black hole temperature for large S-AdS 
black holes. All fits have been obtained on a set of numerical data with black hole radii varying between r + = 10 and r + — 100 
(see text). In parentheses, we show the difference between a given fitting coefficient and the fitting coefficient of the preceding 
mode: as the mode order grows, modes tend to become more and more equally spaced. The axial modes have been labelled 
starting from one to take into account the existence of a purely damped mode. 



Scalar Q = 0) Axial (I = 2) Polar Q = 2) 



Mode 


ujr/T uji/T 


Mode 


ujr/T uji/T 


Mode 


ujr/T uji/T 





7.747 11.157 


1 


7.748 11.157 





7.750 11.153 


1 


13.243 20.592 
(5.495) (9.435) 


2 


13.243 20.592 
(5.495) (9.435) 


1 


13.246 20.585 
(5.496) (9.432) 


2 


18.702 30.021 
(5.458) (9.429) 


3 


18.702 30.021 
(5.459) (9.429) 


2 


18.705 30.010 
(5.459) (9.425) 


3 


24.151 39.447 
(5.450) (9.426) 


4 


24.152 39.447 
(5.450) (9.426) 


3 


24.154 39.434 
(5.449) (9.424) 



with I = (scalar case) and I — 2 (gravitational case). Notice that the labelling of axial QNM frequencies starts 
at n = 1, since we regard the pure imaginary mode found in 0] as the fundamental. All fits have been done on a 
table obtained performing explicit calculations for black holes with horizon radii varying from r + = 10 to r + = 100 
in steps of Ar+ =0.1. Values in parentheses indicate the spacing in the fitting coefficient between a given mode and 
the previous overtone: as the overtone number increase, these numbers tend to a constant. It may also be noted that 
fitting coefficients in this large black hole limit are essentially the same for all three kinds of perturbations. 

In brief, from Table ^ we can deduce a very simple "generalized formula" : for the fundamental mode and the 
overtones of scalar and gravitational S-AdS perturbations, QNM frequencies are well approximated by the fitting 
formula 

lor ~ (7.75 + 5.46n)T, ~ (11.16 + 9.43n)T. (22) 

A relabelling of the integer n is necessary to take into account the existence of a purely imaginary axial gravitational 
mode, but this is just a matter of convention. 

3.1.2. Purely damped modes 

The study of electromagnetic and gravitational quasi- normal modes for S-AdS black holes has unexpectedly revealed 
the existence of axial gravitational modes with pure-imaginary frequency 10]. These modes have the interesting 
property that their frequency does not scale linearly with the black hole radius, but rather with the inverse of the 
black hole radius: an explicit fit (performed using data in the range r + = 10 to r + = 100 at steps Ar + = 0.1) shows 
indeed that for purely damped axial modes with I = 2 (see Table III in 10]) ojj — 1.335/r+. 

In other words, in this case wj scales with the inverse of the mass, as it does for Schwarzschild black holes in flat 
spacetime: hence these modes are particularly long lived. Cardoso and Lemos managed to prove stability for other 
kinds of perturbations, but not for these modes. Indeed, the proof relies on the radial potential being positive, while 
V dd, given by formula (|19fl . can be negative for some values of I and M. 

Electromagnetic perturbations of large S-AdS black holes are also characterized by a set of pure imaginary modes, 
which however scale in the "ordinary" way (u>i ~ r + ). 

We point out that we get essentially perfect agreement with all the numerical results shown in [T^j, with the 
exception of their Table II (lowest QNM of electromagnetic perturbations for 1 = 2). The values we obtain in this 
particular case are given in Table [HI More importantly, we find no evidence for the axial "algebraically special" 
mode with r+ = 1 and pure- imaginary frequency uii = 2, indicated as a "dubious" result in Table III of |10| . In 
both formulations of the eigenvalue problem, our two-dimensional Miiller root-searching routines fail to converge 
there, and even a one-dimensional search on the axis lor = does not show any root. This is not too surprising. 
This algebraically special mode (if it exists) would be characterized by boundary conditions different from those of 
ordinary QNM's: there would be no radiation going down the black hole horizon (To). 

It is well known from the asymptotically flat case [2(j that, as we "turn on" the charge, electromagnetic and 
gravitational modes become "mixed" : as we noted earlier on, only when Q — > the potentials reduce to pure 
gravitational perturbations of Schwarzschild, and the Vy potentials reduce to pure electromagnetic perturbations. 
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TABLE II: Lowest QNM of S-AdS, pure electromagnetic perturbations for 1 — 2. In parentheses we give the values listed in 
Table II of the paper by Cardoso and Lemos. 



r+ uj r (o;g L ) 

0.8 3.224 (2.501) 

1 3.223 (2.496) 

5 3.090 (0.822) 

10 

50 

100 



LP I (CJ? L ) 

0.996 (1.176) 
1.384 (1.579) 
9.822 (10.309) 
16.623 (15.755) 
75.269 (75.139) 
150.133 (150.069) 



TABLE III: Location of minima and maxima of ur(Q) for the fundamental mode. Values listed correspond to some selected 
values of the horizon radius r+ and to different classes of perturbations. A dash ("-") means that our numerical method fails 
to converge before reaching the corresponding minimum or maximum in Q. 

r+ Scalar (I = 0) Axial (I = 2) Polar (I = 2) 



r+ 


mm 


max 


mm 


max 


mm 


max 


100 0.366 


0.474 


0.366 


0.474 


0.366 


0.474 


50 


0.366 


0.474 


0.366 


0.474 


0.366 


0.474 


10 


0.367 


0.475 


0.368 


0.476 


0.369 


0.477 


5 


0.372 


0.480 


0.376 


0.483 


0.376 


0.488 


1 


0.468 


0.571 






0.503 





We will see in the next section that both modes (i.e., those which can be classified as "purely electromagnetic" and 
"pure gravitational" in the S-AdS, Q = limit) show an interesting behaviour as the charge increases. In particular, 
purely damped modes of both classes behave in a very peculiar way: their damping seems to go to infinity (cjj — ► 0) in 
the extremal black hole case, suggesting the possibility that extremally charged RN-AdS black holes may be marginally 
unstable to electromagnetic and gravitational perturbations. 



3.2. Reissner-Nordstrom-anti-de Sitter black holes 



3.2.1. "Ordinary" quasi-normal modes 

A first study of scalar perturbations of RN-AdS black holes was carried out by Wang, Lin and Abdalla ■ The 
authors found that "switching on the charge" results in a breakdown of the linear relations u)r(T), u>i(T) (and in a 
similar breakdown of the linear relation between T and r + ) which is larger the larger is the charge. In a following paper, 
time evolutions of the scalar field wave equation were used to cross-check these results 0] and to study the late-time 
behaviour of the perturbing field. Both numerical methods used in our work (see the Appendix) are frequency-domain 
methods, and break down for large values of the charge, for reasons related to the radius of convergence of power-series 
solutions of the wave equation |12| . However, time evolutions do not suffer of this problem, and were used to show 
that the imaginary part of the frequency attains a maximum at some critical value of the charge Q cr it, to reduce 
again as Q — > Q ex t- This is an important point on which we will comment later. 

Our numerical results for the scalar case agree with those presented in 0, 0] (for example, we are in excellent 
agreement with Table II in ^H), and confirm a breakdown of the linear relations (1221) as the black hole charge 
increases. However, in |l2j it was not only claimed that u>i increases with Q, but also that u>r decreases with Q, so 
that, in their words, "if we perturb a RN-AdS black hole with high charge, the surrounding geometry will not "ring" 
as much and long as that of the black hole with small Q". Our results do not confirm this statement. 

This can be seen, for example, from figure ^ The continuous wiggling "tails" departing from the scalar- mode line 
are numerical solutions of the scalar eigenvalue problem for selected horizon radii and increasing charge. From the 
plot it is clear that, as the charge increases, the imaginary part of the frequency increases, according to the predictions 
in [T^ : however, the real part of the frequency does not show a monotonically decreasing behaviour! 

The situation is clarified in figure|21 where we also consider polar and axial perturbations. Consider first the top two 
panels in figure [3 where we plot the QNM frequency u>r as a function of the normalized charge Q = Q/Q e xt- Each 
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TABLE IV: Fit parameters for the charge dependence of the real part of the QNM frequency for some selected horizon radii 
and different classes of perturbations - see formula 12311 . The accuracy of these fits, which are only valid for Q < 1/3, is 
typically of some parts in a thousand. 

r+ Scalar (I = 0) Axial (j = 2) Polar (I = 2) 
r+ a b a b a b 

100 1.6944 3.6144 1.6943 3.6135 1.6940 3.6107 
50 1.6939 3.6109 1.6935 3.6074 1.6923 3.5963 
10 1.6775 3.5021 1.6771 3.4198 1.6378 3.1729 
5 1.6284 3.1911 1.5882 2.9204 1.4695 2.2189 

TABLE V: Fit parameters for the charge dependence of the imaginary part of the QNM frequency for some selected horizon 
radii and different classes of perturbations - see formula (12 H . The accuracy of these fits, which are only valid for Q < 1/3, is 
typically of some parts in a hundred. 

r+ Scalar {I = 0) Axial (I = 2) Polar (I = 2) 
r+ c d c d c d 

100 0.56827 5.8567 0.56829 5.8556 0.56892 5.8550 
50 0.56839 5.8524 0.56847 5.8483 0.57096 5.8458 
10 0.57193 5.7182 0.57378 5.6222 0.63597 5.5560 
5 0.58145 5.3361 0.58634 5.0124 0.83131 4.7163 

panel corresponds to a different value of r + . As we could expect from the uncharged black hole limit, the scalar and 
mixed gravitational/electromagnetic perturbations are almost isospectral for large black holes, but the isospectrality 
breaks down for intermediate-size holes. Quite interestingly, QNM frequencies show indeed some kind of "damped 
oscillations" as Q increases. In Table Hill we list the location of minima and maxima in lor{Q) for the three different 
kinds of perturbations, and for selected values of r+. For large black holes (which are of greater interest in view of 
the AdS/CFT correspondence), the values of Q for which the frequency attains the local minimum and maximum are 
independent of the black hole radius: they correspond, respectively, to Q = 0.366 and Q = 0.474. 

For small values of the charge, Q < 1/3, an excellent fit of the real part of the frequency at fixed r+ (errors are 
typically of order 10 or smaller, especially for the low-order modes), valid for all kinds of perturbations, is provided 
by a simple polynomial relation: 

lor = u4 0) (1 - aQ 2 - bQ 4 ) , (23) 

where we have denoted by lo^P the real part of the QNM frequency in the zero-charge (S-AdS) limit. The fit parameters 
a and b are given for selected values of r + in Table IIVI Again, they are largely independent of r + and of the kind of 
perturbation considered, at least in the large black hole limit. 

Let us now turn to the bottom two panels in figure |3 For the values of the charge allowed by our numerical 
methods, wj increases monotonically as Q increases. The second derivative of oji(Q) changes sign when the real part 
of the frequency attains the local minimum, and then again at the local maximum . A good fit in the region Q < 1/3 
(errors being of the order of a few percent for the fundamental mode, and of order 1 % or less for the overtones) is 
again obtained using a polynomial relation: 

Lu I =Jp(l + cQ 2 +dQ A ), (24) 

where we denote by tOj the imaginary part of the QNM frequency in the S-AdS limit. The fit parameters c and d 
are listed in Table IVl As for the fitting coefficients of the real part, c and d are largely independent of r+ and of the 
kind of perturbation considered in the large black hole limit. 

Summarizing, for large black holes with charge Q < 1/3, QNM frequencies and dampings are approximated within 
about 1 % by a very simple formula: 



u R = (7.75 + 5.46n)T (0) (l - 1.694Q 2 - 3.61Q 4 ) , (25) 
uji = (11.16 + 9.43n)T (0) (l + 0.568Q 2 + 5.86Q 4 ) , 
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where — (1 + 3r^_)/47rr + is the S-AdS (uncharged) black hole temperature. Notice that, for fixed Q, the 

dependence of lur and ujj on the angular index I is extremely weak, so this formula holds with good accuracy also for 
angular multipoles larger than the fundamental. 

3.2.2. Purely damped modes: are extreme RN-AdS black holes marginally unstable? 

Modes which become "purely electromagnetic" or purely imaginary in the S-AdS limit deserve a separate discussion. 

Let us first consider the pure-imaginary modes for axial gravitational perturbations. In figure |3] we "track" the 
imaginary part of these modes, starting from the corresponding zero-charge limit, for different values of the horizon 
radius. As explained when we discussed purely imaginary QNM's for S-AdS black holes, our numerical results do not 
confirm the existence of an algebraically special mode in the zero-charge limit at r+ = 1: as we decrease r+, the root 
finder breaks down when r + ~ 1.04 (dashed line in the plot), corresponding to a S-AdS frequency uji = 1.860. 

It is clear from the plot that all modes show a tendency to approach the horizontal axis as the charge approaches 
the extremal value. This tendency is much clearer if we use a linear scale (instead of the log scale used in this plot, 
which is necessary to display all modes together). Unfortunately, for the convergence reasons discussed earlier, we 
cannot push our calculations to the extremal limit. However, our plots lead us to conjecture that all pure-imaginary 
axial modes are such that wj — > as Q — > Q ex t ■ 

Our conjecture finds firmer ground in the behaviour of the modes for the V{~ potential which reduce to purely- 
damped electromagnetic perturbations with I = 1 and / = 2 in the S-AdS limit. In this case, we were able to track 
roots for all values of the charge Q < 1, and our results are shown in the two panels of figure 0] the imaginary part 
of purely-electromagnetic modes in the Schwarzschild limit does indeed tend to zero as Q — > 1. 

For small black holes, such a behaviour is not observed. Small black holes in S-AdS do not have pure imaginary 
electromagnetic modes. Tracking again the modes of the V{~ potential (figure [SJ we see that their real part decreases 
as a function of Q, but the imaginary part does not: indeed, u>i(Q) is decreasing when r + > 2, but it shows a relative 
minimum when r + = 1.5, and is roughly constant or slowly increasing when r + < l! 

Our results suggest the possibility that, unless the mode amplitude tends to zero in the extremal limit, large extremal 
RN-AdS black holes are marginally unstable to electromagnetic (and perhaps, axial gravitational) perturbations. This 
result may be extremely interesting in view of the AdS/CFT conjecture. 

We only make a couple of comments at this stage. First, notice that even for S-AdS black holes, a proof of stability 
to large black hole, axial gravitational perturbations does not exist [ljj. Second, as we recalled earlier, time evolutions 
of scalar fields in RN-AdS spacetimes have shown that, for Q greater than some critical value Q cr it, is a decreasing 
function of charge |15| . The authors suggested a possible connection between this behaviour and the existence of a 
second-order phase transition for extremal RN-AdS black holes |2j| • Development of more efficient numerical methods 
is needed to verify such a connection between QNM's and thermodynamical phase transitions of extremal RN-AdS 
black holes. An understanding of the link between our study and the dynamical instability discussed by Gubser and 
Mitra |24l |25| , which in their words persists "for finite size black holes in AdS, down to horizon radii on the order 
of the AdS radius", could shed further light on the AdS/CFT conjecture. We notice that an interesting discussion of 
the connection between dynamical and thermodynamical stability in AdS spacetimes can be found in [26j. 

4. CONCLUSIONS 

We have studied scalar, electromagnetic and gravitational perturbations of a RN-AdS spacetime, and computed its 
QNM's. Our results extend previous computations carried out for scalar, electromagnetic and gravitational pertur- 
bations of S-AdS black holes and for scalar perturbations of RN-AdS black holes For "large" black holes, 
whose horizon is much larger than the AdS radius, different kinds of perturbations are almost exactly isospectral; this 
isospectrality is broken when the black hole's horizon radius is comparable to the AdS radius. Contrary to previous 
claims, we find that, as the black hole charge Q increases, the real part of the QNM frequency does not continuously 
decrease, but rather shows a minimum followed by a maximum. The location of these extrema as a function of 
Q = Q I Qext is independent of the black hole radius for large black holes: minima are located at Q — 0.366, and 
maxima at Q = 0.474. We found that the imaginary part of the QNM frequency Wj is generally, at least for small 
values of Q, an increasing function of the black hole charge, whose second derivative changes sign corresponding to 
the minima and maxima of u>r(Q). We also obtained the very accurate fitting formulas (|25|l . which are valid for large 
black holes having charge Q < Q ex t/3- 

A very interesting result concerns electromagnetic and axial perturbations of large black holes. These perturbations 
are characterized by the existence of purely damped modes. We advanced numerical evidence that the damping of 
these modes approaches infinity as the black hole charge approaches the extremal value. If the corresponding mode 
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amplitude does not tend to zero, this would imply that extremally charged RN-AdS black holes are marginally unstable. 
Our results are relevant in view of the AdS/CFT conjecture, since, according to it, the AdS QNM's yield thermalization 
timescales in the corresponding conformal field theory. It would be desirable to relate our results to the existence of 
a second-order phase transition for extremal RN-AdS black holes and to a dynamical instability which has been 
studied in the past [2J, |25j , and which also holds for finite size RN-AdS black holes having horizon radii larger than 
the AdS radius. We believe that investigating these connections could clarify many issues related to the AdS/CFT 
correspondence. 
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APPENDIX A: NUMERICAL COMPUTATION OF QUASINORMAL MODES 
1. The Horowitz— Hubeny method 

Here we briefly describe the first numerical method we used to compute the quasinormal modes. This method has 
been described elsewhere 

HQHH3, so we only give the essential formulas, which encompass the slightly more general 
case considered here. For a generic wavefunction obeying one of the five wave equations JSJ), l|lUfl or (|15|) . write 

x (,(r) = e 4ur *Z{r) (Al) 

and get the equation for ip(r) 

f( r )^T + t/'W " 2 H ^ " V(r)i>(r) = 0, (A2) 

where V(r) = V(r)/ f(r). Introducing a new variable x = 1/r to "compactify" the region outside the black hole, and 
defining x + = l/r + , this equation is rewritten as 

s(x)^ + t{x) + T^K^) = (A3) 

ax A x — x + ax [x — x+) 



where 



r x 5 - x 4 - x 2 - Q 2 x e 

s(x) = , (A4) 

x — x+ 

t{x) = 3r a; 4 - 2a; 3 - 4QV - 2iux 2 , (A5) 
u{x) = (x-x + )V(x), (A6) 



and 



X + 

Now look for a series solution in the form 

oo 

4>{x) = ^a k {x-x + ) k . (A8) 

k=0 

This series solution obeys the boundary conditions required for quasinormal modes if the condition 

oo 

7/K0) = 5> fe (-x+) fc = (A9) 

k=0 
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is satisfied. The coefficients a k = afc(w) can be obtained through the following recursion relation: 

1 £1 

cij = -— ^ - ft + ktj-k + Uj-k] a k . (A10) 

In this formula the s^'s, i^'s and u k s are coefficients in a Taylor series expansion around x = x+ of the functions 
s(x), t(x) and u(x) respectively, and 

P j =j(j-l)s +jt . (All) 



2. The Frobenius method 

As a check of our results, we have applied an alternative method, originally suggested by Leaver [Tsj and applied in 
[ll| to the axial perturbations of S-AdS BHs. We have generalized the approach to RN-AdS BHs. Since this method 
is computationally much faster, and very easy to implement for the axial potentials V[~ and V 2 ~, we give the details 
here. 

Rewrite equation (|A3(1 multiplying it by (x — x+)/x 2 to find: 

, .d 2 ib(x) , ,dib(x) V(x) , , , . . „. 

p^^k 1 + + ~^^ x) = (A12) 

where as usual V = V/ f , and furthermore 

p(x) = {2Mx 3 - x 2 - 1 - QV), 

q(x) = (6Mx 2 -2x~4Q 2 x 3 -2iu). (A13) 
Now look for the solution in the form of a Frobenius series which is regular at the event horizon: 

u(x) = J2"n(^ ± ) , (A14) 



where ciq = 1. 

In the RN-AdS case, as is the case for Reissner-Nordstrom black holes in asymptotically flat spacetimes the 
expansion coefficients are determined by a four-term recurrence relation of the form: 

a Q ai + (3 a = 0, 

a\a 2 + Piai + 71 a = 0, 

a n a n+ i + /3„a„ + 7„a„_i + 5 n a n -2 = 0, re = 2,3,... (A15) 

where the recursion coefficients are given in terms of the parameters appearing in the black hole potentials V~ , defined 
in equation (|T3)) . as follows: 

a n = 2(n + 1) [{n + 1)(1 - 3Mx + + 2Q 2 x\) + ioj/x+] , 

p n = n{n + 1) (6Mx+ - 1 - 6Q 2 x 2 + ) + 1(1 + 1) - p 3 x+ - ^-x\, 

2„,2 o \ , „ „ . PiPj 2 



7 „ = (n 2 - 1) (4Q 2 a^ - 2Mx+) + Pj 



S n = -(n-2)(n + l)Q 2 x 2 + -^x 2 + . (A16) 

In the zero-charge limit Q = we recover a standard, three term recursion relation [Tl|. Once the a„'s are known 
(in practice, up to some finite value of n determined by the accuracy we want to achieve), the QNM's are determined 
by imposing the condition u(x = 0) = 0. 
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FIG. 1: Plot of the S-AdS QNM frequencies for the fundamental scalar mode with I = (continuous line), for the first 
non-purely damped axial mode with / = 2 (dashed line) and for the fundamental polar mode for I = 2 (dotted line). All 
calculations are started at r+ = 100 (top right of the diagram). The step in r+ is reduced by Ar+ =0.1 and the calculation 
is iterated until our numerical code fails to converge (bottom left). Corresponding to selected values of the horizon radius 
(namely, r+ = 100, 50, 10, 5, 1, 0.8, 0.6, 0.4, whose S-AdS frequencies are tabulated in Table I of Horowitz and Hubeny), we 
"switch on" the charge and follow the modes in the complex plane. The resulting trajectories are the "tails" departing from 
the continuous scalar QNM line. Axial and polar perturbations are almost isospectral both in the large and in the small black 
hole limit. 
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FIG. 2: The top two panel shows the characteristic "wiggling" of the real part of the QNM frequency, ujr, as a function of the 
normalized charge Q = Q/Q e xt, for selected values of the horizon radius r + ; the bottom two panels show the imaginary part 
of the QNM frequency, u)i(Q), for the same black holes. Scalar perturbations are denoted by solid lines, axial perturbations 
by dashed lines, and polar perturbations by dotted lines. The isospectrality of different kinds of perturbations, which holds in 
the large black hole limit, is clearly lost as the black hole "size" becomes comparable to the AdS radius, r + ~ 1. Notice that 
W "(Q) = when u>' R (Q) = (a prime denoting differentiation with respect to Q). 
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FIG. 3: Imaginary part of the purely damped mode reducing to pure axial gravitational perturbations with I — 2 in the 
zero-charge limit. Starting from S-AdS results, we track the modes for selected values of the horizon radius r+ (indicated to 
the right of each curve). The calculation is terminated when our root finder fails to converge. The dashed line indicates the 
last mode we can find before the numerical method breaks down as we approach the "algebraically special" (wj = 2, r+ = 1) 
frequency. 
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FIG. 4: Imaginary part of the purely damped mode reducing to pure electromagnetic perturbations with I = 1 (left panel) 
and / = 2 (right panel) in the zero-charge limit. Continuous, dashed and dotted lines refer, respectively, to black holes having 
horizon radius r+ = 100, 50 and 10 (top to bottom in the plots). For / = 1, we also plot the pure imaginary mode with r+ = 5 
(dash-dotted line). 




FIG. 5: Real (left panel) and imaginary parts (right panel) of the small-black hole QNM frequencies for modes reducing to 
pure electromagnetic perturbations with / = 1 in the zero-charge limit. The size of the black hole horizon radius corresponding 
to each curve is indicated in the inlay. 



